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Abstract. After a successful development of theoretical and numerical works on Fermi acceleration at relativistic shocks, 
some difficulties recently raised with the scattering issue, a crucial aspect of the process. Most pioneering works were 
developed assuming the scattering off magnetic fluctuations as given. Even in that case, when a mean field is considered, 
its orientation is mostly perpendicular to the shock normal in the front frame, and this tends to quench the scattering 
process. Solving this difficulty leads to address the issue of the generation of very intense magnetic fluctuations at short 
wave lengths. The relativistic motion of the shock front let the cosmic rays to visit upstream during a very short time 
only, making this generation of magnetic fluctuations very challenging. Anyway there is some hope to solve the problem. 
Thanks to a recent work by Spitkovsky (2008) |18], we know that the process works without any mean field and now we 
have to investigate up to which intensity the mean field can be amplified for allowing Fermi process with appropriate fast 
instabilities. In this presentation, the coUisionless shock structure in relativistic regime is sketched, the scattering issue is 
presented, and the instabilities that can provide the expected magnetic field amplification are presented as well. Although 
there exists observational evidence that particles are accelerated in relativistic flows and are distributed according to a power 
law suggesting a Fermi process, the drastic conditions for Fermi process to work are not always clearly fulfilled. 
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INTRODUCTION 

Important progress were accomplished during this last 
decade about the acceleration of particles governed by 
Fermi process at relativistic shocks 01], H], lH, lH]. 
In these works a prescribed scattering process was in- 
troduced in the calculations and in the numerical works. 
In Lemoine & Pelletier 03 [12] a scattering medium 
was simulated by introducing magnetic fluctuations with 
Fourier modes distributed according to a power law 
spectrum. The same results were obtained, in particular 
the spectrum index between 2.2 and 2.3; however with 
the assumption that the pitch angle only is involved in 
the scattering process; which was also assumed by the 
previous works mentioned. More detailed simulations 
by Niemiec, Ostrowsky |10], then Niemiec, Ostrowski 
and Pohl lHHl . and also by Lemoine, Revenu llTill . then 
Lemoine, Pelletier, Revenu lll3ll revealed the difficulties. 

• The mean field suffered by particles is essentially 
transverse to the flow in the front frame, if the angle 
Qb of the mean field with respect to shock normal 
be larger than 1 /Ys, which is obviously the most 
frequent case. 

• The downstream magnetic field is very compressed 
and weakens particle diffusion, therefore the univer- 
sal spectrum is not obtained [10] [lUl 13] • 

• An usual turbulence cascade spectrum from large 
scales to small scales cannot provide the expected 



scattering d [Hsl] , 

Finding circumstances favorable for Fermi process at a 
relativistic shock is very challenging. 

Because it is useful for these issues, we first indi- 
cate what a coUisionless shock structure is, especially in 
proton-electron plasma with a quasi perpendicular mean 
field. Then we present the scattering issue and the strong 
challenge for having Fermi cycles. Then we propose sev- 
eral fast instabilities that are good candidates for am- 
plifying the magnetic fluctuations that would allow the 
expected scattering. However we found stringent condi- 
tions on the mean field for that amplification. 

The recent remarkable work by A. Spitkovsky [ lisl] 
shows through PIC simulations, that the Fermi process 
works at relativistic shocks in a non-magnetized pair 
plasma. The shock structure is quite different. The mag- 
netic fluctuations are produced by a Weibel instability. 
We will indicate how weak a mean field must be for the 
instability to be operative even with a proton-electron 
plasma. We propose also another instability that works 
with a larger mean field, however still quite weak. 

There is observational evidence in Gamma Ray Bursts, 
especially for the afterglows, that both the magnetic field 
has been amplified upstream and that power law distribu- 
tions of high energy particles has been formed suggest- 
ing a Fermi process (see, for instance Li & Waxman 06 
lUTT ). Moreover GRBs often explode in Wolf-Frayet star 
winds, which are much more magnetized than the inter- 



stellar medium (Ud-Doula et al. 08). 

PHASE SPACE LOCKING BY THE 
TRANSVERSE MEAN FIELD 

How a coUisionless shock is set up? 

A coUisionless shock is built with the reflection of a 
fraction of incoming particles at some barrier, generally 
of electrostatic or magnetic nature, except in the case of 
a pure electron-positron plasma without any mean mag- 
netic field. It suffices that some downstream obstacle 
constrains the flow to slow down and thus to compress. 
In a proton-electron plasma carrying an oblique magnetic 
field a barrier of both electrostatic and magnetic nature 
raises. Because the magnetic field is frozen in most part 
of the plasma, its transverse component is amphfied by 
the velocity decrease; thus a magnetic barrier is set up. A 
fraction of the incoming protons is reflected back. A po- 
tential barrier raises also because the electron distribution 
is close to Boltzman equiUbrium and thus ~ Tglog ^ 
with an electron temperature that grows to a value com- 
parable, but likely different, to that of protons, which 
reaches Tp ~ (Fj — \)mpC^. The potential barrier grows 
up to a value that allows reflection of a significant part of 
the incoming protons, which means eA<t> ^ (Fy — l)mpC^ 
(we assume an ionic population mostly composed of pro- 
tons). This potential barrier reflects also a fraction of pro- 
tons, but favors the transmission of electrons that would 
otherwise be reflected by the magnetic barrier also. The 
reflection of a fraction of the protons insures the matter 
flux preservation against the mass density increase down- 
stream. However because the magnetic field is almost 
transverse, an intense electric field E = jS^fi energizes 
these reflected protons such that they eventually cross 
the barrier. Interactions between the different streams of 
protons generates a turbulent heating of the proton pop- 
ulation, which takes place mostly in the "foot" region. 
The foot region extends from the barrier upstream over a 
length If = rL,F = r^V'^/to^^F. measured in front frame, 
and measured in upstream frame, it leads to 



tF,u 



COcirs sin Ob 



(1) 



We assume that the field is almost perpendicular in the 
front frame, but take into account of an angle Ob in the 
co-moving upstream frame such that sin 9b > I /F^. The 
downstream flow results from the mixing of the flow of 
first crossing ions (adiabatically slowed down) with the 
flow of transmitted ions after reflection. All the ingredi- 
ents of a shock are reaUzed. The three ion beams in the 
foot interact through the "modified two stream instabil- 
ity"; this constitutes the main thermalisation process of 



the ion population that continues with only two beams 
passed the potential barrier. These anomalous heating 
processes require an appropriate kinetic description. 

Ramp, overshoot, foot. 

Entropy production in the shock transition comes from 
two independent anomalous (caused by coUisionless ef- 
fects) heating processes for electrons and ions. In non- 
relativistic shocks, electrons reach a temperature larger 
than ions; however we do not know yet whether this is 
still the case in relativistic shocks. Moreover they ex- 
perience a supplementary heating with the convection 
electric field. Whereas the growth of ion temperature de- 
velops on scale £f, electron temperature grows on the 
very short scale £,■; this defines the "ramp" of the shock. 
Actually, because of this strong gradient where an in- 
tense transverse electric current is concentrated, anoma- 
lous heat transfer occurs through the ramp. This electron 
heating is described by Ohm's law in the direction of the 
convection electric field (z-direction), namely 



ricdB 
An dx 



with£ = j3,5o ; 



(2) 



it indicates that the variation scale is the relativistic re- 
sistive length given "oy ir = ^ = 8e^-^ {^x < 0). This 
is a very short scale comparable to the electron inertial 
length 8e when the anomalous resistivity is so strong that 
the effective collision frequency is of order (Ope. This is 
the scale at which the electron population experiences 
Joule heating, this is thus the growth scale of three major 
quantities, namely, the potential, the magnetic field and 
the electron temperature. Even if the scale 5e is estimated 
with electrons so hot that they have the relativistic mass 

TsTHp (i.e. 5e^ 8 = ( ^^'"^^2 y^^)-, it remains smaller than 
the foot length; indeed 
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^ < 1 <^ < npTsmpC^ ■■ 



(3) 



that strong inequality is a natural requirement for a rel- 
ativistic strong shock. By the way, the condition for the 

Fermi process to be operative will turn out to be more 
restrictive for the field intensity, as will be seen further 
on. 

This strong Joule heating results from the very in- 
tense cross-field current and, actually, is produced by 
microturbulence driven by the current itself. Probably an 
anomalous diffusion of electron temperature occurs that 
smoothes the temperature profile; however it is not yet 
identified in relativistic shocks. 

The field profile can be obtained by prescribing a 
velocity profile decreasing from 1 to 1/3 over a distance 



much larger than £r. The profile displays a ramp at 
scale followed by an "overshoot" before reaching the 
asymptotic value 3Z?o- Electrons have to jump this barrier 
in order to come from downstream to upstream. The 
formation of a foot over a length of order £f much larger 
than the resistive length, which characterizes the ramp 
and the overshoot, is the main feature of collisionless 
shocks, perpendicular or oblique (see figure(l).) 

In the case of an electron-positron plasma, when 
a magnetic field is considered, no electrostatic barrier 
raises, only the magnetic barrier occurs. When the mean 
magnetic field is very negligible, a barrier can raise only 
through the excitation of waves. This is a crucial as- 
pect of the simulation of A. Spitkovsky |18] that shows 
that the Weibel instability can produce particle reflec- 
tion through the wave growth that exerts a ponderomo- 
tive force. 

When electron reach energies comparable to FsinpC^, 
the situation becomes similar to the electron-positron 
plasma shock. 

The strucutre is described with two scales 5 and If 
and three small pai-ameters: ^cr, the fraction of incoming 
energy density converted into cosmic ray pressure at 
front, £,B the ratio of magnetic energy density over the 
incoming energy density and j!- . 



Particle motions 

Particle motions upstream 

Besides the part of "thermal" protons tiiat are reflected 
back on the electrostatic barrier and eventually cross 
the ramp after having gain an energy of order r^mpC^, 
some amount of the downstream protons comes back 
upstream, scatters off the magnetic field and comes back 
downstream with an energy of order F^mpC^. This is 
the first cycle of Fermi acceleration, which is always 
possible for a fraction of the particles and that generates 
a primary population of high energy cosmic rays. The 
question is to know whether further Fermi cycles are 
possible. 

Let us first look at the kinematics of cosmic rays up- 
stream in a large scale mean field B = Bey parallel to 
the shock plane; it is convenient to make this investiga- 
tion in the co-moving frame, where the electric field van- 
ishes and where particle orbits are circles characterized 
by the Larmor radius. The momentum py is conserved. 
With simple geometrical arguments one can derive the 
gyro-phase variation as a function of the initial angle 
00 of the momentum in the vertical plane with respect to 
the flow (see figure 2). 



One finds 



sin(0() + V^) - sin Oq = p^W ■ 



(4) 



The angle is necessarily such that cos Oq > jSj and thus 
I sin 00 1 < 1 /^s- The small angle approximation leads to 



1 



1 



-V/" + 0or+0o = ^ 



(5) 



The maximum deviation is obtained for 9q — and 
equals ± The protons are deviated upwards, whereas 
electrons are deviated downwards. One easily checks 
that no deviation occurs at the maximum entrance angle 
±1 /Fs. The visit time is thus given hy y/ = (OL^t ^ 1 /F, 
and thus Af ^ ti/Ts- The penetration length of the cosmic 
rays into the upstream flow is therefore of order ri/Ys- 
The distance between the orbit and the shock front is 
given by Ax = '"t (sin ( 0o + — sin Qq — j3s i/) and reaches 
a maximum for sin(i^ + 0o) = f- and thus Ax„,a.v ^ p-; 
this indicates the scale at which magnetic perturbations 
could scatter the cosmic rays. 

This relation allows to determine the energy gain after 
a Fermi cycle dud as a function of the entrance angle. 
Using back and forth Lorentz transformations we obtain 
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where we inserted the relative velocity between the up- 
stream and the downstream flows (actually in the lat- 
ter formula we inserted a converging velocity — j3r and 
counted j3r positively): 



j3, = (A,-j3^)/(l-j3„j3,/) 



1 - 



1 

F? 



The angle 0o = produces the maximum gain: G(0) = | 
and the minimum gain, G = 1, is obtained for 0o = 
± 1 /F.V. The gain monotonically decreases between these 
two values. The gain averaged over this angle interval is 
close to the value 2. The large gain of order F^ cannot 
be reached with particles coming from downstream to 
upstream; only the first Fermi cycle udu can produce 
such an extreme gain (see 101, iQ]). 



Particle motion downstream, the "no return " sector 

A similar dynamical system governs downstream mo- 
tion; j3s is replaced by /3(/ ~ 1/3 for a strong relativistic 
shock. An incoming proton with jS^ > and jS^ < ^id has 
no possibility to come back to the shock. Only intense 
magnetic perturbations on short scales (on scale com- 
parable to the Larmor radii of the particles measured in 
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FIGURE 1. Shock structure: This drawing roughly sketch the profile of the main physical parameters that varies in the foot, the 
ramp, the overshoot. 
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FIGURE 2. Particle motion upstream 



front shock) could produce efficient scattering that would 
allow further Fermi cycles. For instance if strong turbu- 
lence makes the magnetic field to reverse on short scales 
this would probably make Fermi cycles possible. 

A simple geometrical argument can explain the "no- 
return" condition (see figure 2). A similar geometrical 
construction can be done to look at the orbit of an incom- 
ing particle that rotates in the co-moving downstream 
flow whereas the shock is moving away at the speed c/3. 

One finds the relation 



sin(0o + W)- sin 0o = ^ or p^W 



(7) 



with the condition that incoming particles are such that 
cos 00 > 1/3. This equation has roots in some interval 
and no roots in the complementary interval. It can be seen 
first by considering two simple cases. First consider the 
case 00 = 0; the equation has no root because > ~ 5 
for all y/. Now consider the case 0o = f ; the graph shows 
evidence of a root. 

A particle crosses the shock from downstream to 
upstream if it has an angle 0^ G (— 0*,0*), such that 
COS0* = 1/3 (0* ~ 1.231 rac/), (a "pitch" angle corre- 
sponds to an angle 0: — ^ on the Larmor circle) which 
corresponds to a = — ;r/2 between —0* — njl and 
0* — 7r/2 on the Larmor circle (a* = —0.349 rat/). Only 
some interval of initial angles allows the particle to re- 
cross the shock d ^ u. This interval starts at 0* and ends 
at Qcr such that Qcr + W ~ ^* ^ this angle is thus the root 
(other than 0*) of the equation 



sin — sin 0* = 



— 0* 



(8) 



The root is between n/2 and n, with 0„- — 1. 7 10 rad. 
The interval (0*, Ocr) defines a sector of no-return. The 
particles that make the cycle dud come back downstream 
with an angle shghtly larger than 0, and thus flow in the 
sector of no-return. It corresponds to angle a^r between 
and n/2 on the Larmor circle, acr = OAAOrad. 



Requirement on turbulence for Fermi cycles 

The previous paragraph was devoted to motions in a 
regular transverse field. However the behavior we have 
analyzed with an ordered large scale field also applies 
for ordinary turbulence, that cascades from large scales 
to short scales. This is an important point that has been 
argued in Lemoine, Pelletier & Revenu [13]. In particu- 
lar the figures in this paper clearly show that the Poincare 
velocity mapping, in a plane parallel to the shock plane, 
is similar between regular motions and irregular motions 
in a Kolmogorov turbulence field and reveals the sector 
of "no return". This explains also the results by Niemiec 



et al. ( ifToll . ifTTll ). In order that particles undergo several 
Fermi cycles and form a power law distribution, intense 
scattering should develop in this flow. Magnetic turbu- 
lence should be excited for the Fermi process to be op- 
erative, upstream at time scale shorter than fi/F^ over a 
precursor of length r^/Fj and downstream at time scale 
shorter than t^ over scale shorter than r^, with 5B at that 
scale more intense than the large scale field... 

Thus it is particularly important to look for the gen- 
eration of intense short scale turbulence upstream, and 
preferentially incompressible turbulence that can reach 
higher level than compressible one. The requirement 
on the relative amplitude A of the magnetic fluctuation 
can be explored by looking at the variation of the an- 
gles and defined by )3vCos0, jS^ — sin cos and 

~ sin sin . 



< A0^ > 



3A^^T..A. 
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(9) 
(10) 



Careful analysis 111511 shows that the amplitude must be 
quite high, namely A > ri/Fsic, £c being the coherence 
length of the magnetic fluctuations, where the Larmor 
radius is estimated in the quadratic averaged field and 
Fermi acceleration is restricted to hmited range of Lar- 
mor radii; 

r,4<rL<Ar,4. (H) 



INSTABILITIES WITH A MEAN FIELD 

The expected instability that would scatter the cosmic 
rays upstream is supposed to do it during the precursor 
crossing. Two states of the precursor must be considered: 
one (case A) is the stationary situation where Fermi pro- 
cess is supposed to be developed such that the precursor 
length is determined by the diffusion length of the most 
energetic Cosmic Rays; the other (case B) is the starting 
situation where the cosmic rays of the first Fermi cycle 
are scattered by the mean field. 

In case A, when particles are scattered off short scale, 
but intense, magnetic fluctuations, the scattering fre- 
quency is 



V, - ce^ < 5B^ > 4/£ 



(12) 



and the diffusion length ^ c/Vs (the diffusion coeffi- 
cient is D = c^/3Vs °^ E^); but the scattering must de- 
velop over a length of order c/ViF, because of the front 
motion. For a temporal growth rate yinsi, the growth fac- 
tor is 

„ YinstC 



(13) 



and the energy of the magnetic fluctuations is amplified 
by a factor ^„.rt ^ e^^""" in the linear regime. The insta- 
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FIGURE 3. Particle motion downstream 



bilities are efficient if 

where Va is the Alfven velocity, = Va/Cci is the 
minimum scale for MHD description. Even MHD insta- 
bilities could develop. 

In case B, (see eq.Q) the growth factor 



Ginst — 



Ti n .S t ^ 



(15) 



a)„ysriSin0B 

and thus the temporal growth rate must be larger than 

}^„,r > «„i3srvSin0B . (16) 

Note that we keep sin 0b > 1 throughout the further 
discussions. This last constraint is the most severe. Any 
MHD mode, whose temporal growth rate is necessarily 
smaller than (Od, cannot be amplified unless the magnetic 
field would be parallel to the flow. Nor any extraordinary 
ionic mode that would have a frequency close to o);,, 
such a mode with a frequency close to the lower hybrid 
frequency would be better. Actually if the scattering 
process develops through any fast enough instability, it 
enlarges the precursor and makes the instability growth 
easier. 

Of course the electronic modes are fast enough, but 
their capability to scatter high energy cosmic rays is very 
unlikely. 



At MHD scales 

For the frequently valid condition ^a^s sin 0b ^ 1, the 
precursor has a length much larger than the minimum 



scale for MHD description {tmhd/^F,u = j3Ar.sSin0B). 
The growth of MHD modes is expected mostly in case 
A. There is an interesting possibility of Alfven mode am- 
plification due to the quasi Tcherenkov resonance that 
selects Alfven modes such that ©a = A;||Va = kxVs- The 
growth rate is 



(P^coLcoa) 



1/3 



(17) 



where cOpt. = (47rncTe^/7*'Wp)'^^ with y^, ^ (there is a 
factor F, because the energy of the particles in the front 
frame is of order FsinpC^, and another F^ because of the 
frame motion). The cosmic ray density tier is measured 
in front frame and is related to the cosmic ray pressure 



'2^2 
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at front frame P„ by P„- = UcrT^mpC^ = ^crPu^c- 
that (^(,- = ncr/niiFs- The frequency (Oa is smaller than 
COci but can be close, especially for the right mode. 

When the magnetic field is almost parallel, i.e. when 
0B < 1 /F,, Bell's instabiUty H] lUl] ||^ (see ^ in rel- 
ativistic regime) can develop. It is excited by the charge 
current carried by the cosmic rays in the precursor. Non- 
resonant waves of wavelength shorter than the Larmor 
radii are exited without any response of the cosmic ray 
plasma. Its growth rate in upstream frame is of order 
^cVa with kc = 4nJcr/cBu- Unfortunately this field con- 
figuration is not generic for relativistic shock. However 
for a quasi perpendicular field, in the same spirit as 
Bell's instability, the electric charge carried by the cos- 
mic rays can excite non-resonant compressive modes for 
which the cosmic rays have no response because the 
waves are at short scales compared to the Larmor radii 
also. The cosmic-ray electric charge is almost neutral- 
ized by the ambient plasma electric charge in the front 
frame and a typical wave number is associated to it: 



A;* = ricreBo/ PuT^c^ (defined in front frame). The spa- 
tial growth rate measured in front frame is (see Pelletier, 
Lemoine, Marcowith 08 11511 ) 

7x - {k%YlHorky<^{hklYI^ (18) 
- (;t*fc,)'''^for;t,> (fc*fc2)i/3 

(x is along the flow and y along the mean field). 

In case A, the growth requirement is likely fulfilled 
since 

rirr 1 



^2 



- , (20) 
n« PA^s Pa 

and is larger than 1 for jS^r^ sufficiently small. In case B 
this is more difficult because k^ip = ^-p- = which 
is rather small, but this can be compensated by the wave 
numbers that can be much larger than k^,. 

The development of the compressive instability in case 
A is not expected to solve the scattering issue. It produces 
more likely heating (see further on for a more detailed 
discussion) and the magnetic field fluctuations remain 
moderate because of the limitation of the density deple- 
tion. However a nonlinear investigation, through numer- 
ical simulations, should be useful to look more deeply in 
the consequences of the instability. 



At micro-scales 

When the ambient magnetic field is disregarded, the 
reflected particles and the fraction of particles that par- 
ticipate to the first Fermi cycle, constitute a relativistic 
cold beam that pervades the ambient plasma and trig- 
ger micro-instabilities. One is the two stream instability, 
which mostly ampUfies the electrostatic Langmuir field 
through a resonant interaction such that co — kxVs = 
(Landau-Tcherenkov resonance) and also an electromag- 
netic component because of the density inhomogeneity. 
The other is the Weibel instability, with k^ = (it is non- 
resonant or the resonance is somehow reduced to O) = 0), 
which is mostly electromagnetic with a low phase ve- 
locity so that the magnetic component of the wave is 
dominant. That instability is suitable for developing par- 
ticle scattering. Its growth rate is of order Ps(Op*, where 
a)p* = {Anncre^/YtinpY^^ with y* ^ T^- 

When the ambient magnetic field is considered, the 
relativistic stream of particles that pervades the upstream 
plasma is also seen as a cold beam from the upstream 
comoving frame, with a transverse dispersion of order 
1 /Fs, but with a finite penetration length (Ifm)- The 
beam is almost non-magnetized and can undergo quasi 
Tcherenkov resonance, as previously, but with waves of 
the ambient plasma where the magnetic field plays a role. 
Depending on the nature of the excited waves, the ambi- 
ent plasma is more or less magnetized with respect to 



the wave dynamics. Ambient electrons are easily magne- 
tized, however it turns out that even magnetized, Weibel 
instability can develop provided that ambient protons are 
not (see [16]). Between the inertial scale, 8e = c/cOpg, of 
electrons and the inertial scale 5,- = c/cOp, of the ions, 
there is an intermediate dynamical range between MHD 
(at scales larger than 5,) and electron dynamics (i.e. dy- 
namics for which electron inertia is relevant), where fast 
enough waves can be excited by the relativistic stream. 
In the continuity of right Alfven waves (the left ones are 
completely absorbed at the ion-cyclotron resonance), for 
quasi parallel (with respect to the mean field) propaga- 
tion there are whistler waves, that are electromagnetic 
waves with a dominant magnetic component. For quasi 
perpendicular propagation, there are the ionic extraordi- 
nary modes, which have frequencies between the ion- 
cyclotron frequency and the low-hybrid frequency (ob- 
tained for large refraction index) and which are mostly 
electrostatic with a weaker electromagnetic component. 

For scattering purpose, the whistler waves are the most 
interesting in this intermediate range, (they are excited in 
the foot of collisionless shocks in space plasmas), and for 
pre-heating purpose, the extraordinary ionic modes are 
more interesting (they are actually used for additional 
heating in Tokamak), as in collisionless non-relativistic 
shocks in solar wind. 



Whistler modes instability 

In usual conditions, whistler waves, that are right 
modes, resonate at Landau-synchrotron resonance with 
electrons only. In this peculiar situation, they resonate 
at the quasi Tcherenkov resonance, as mentioned above 
(see ||T6] for details on these developments). Therefore 
all relativistic particles participate in the interaction and 
thus experience scattering if the waves are excited by the 
stream. Actually a fast instability results from the res- 
onance with the following characteristics: a frequency 
(On = (Ocik^df = (Ocek^dg, kx = COvv/c < k and a growth 
rate 

Yin,,^{fi}(£>%(0,,)"\ (21) 

within some unimportant angular factors. This instability 
exists in presence of a mean field and its growth rate 
looks a little faster than the Weibel growth rate when 



An^"^'^ ^ F? 



(22) 



However, when the field satisfies this inequality, the 
Weibel instability is even quenched (see further on) and 
the instability of whistler modes alone will grow both in 
case A and B. 



Instability of extraordinary ionic modes 



At large wave-lengths, the extraordinary ionic waves 
become the magneto-sonic waves of MHD description. 
They are excited by the resonant interaction with the 
beam. The growth rate is a sizable fraction of the mode 
frequency cOx, with cOd < COx< COu, (where the low-hybrid 
frequency is such that (ofi^ ~ O^,- + (0^^, with the assump- 
tion (Opi ^ (Oci). Whereas a non-relativistic beam selects 
a resonant mode of high refractive index, which corre- 
sponds to a high frequency close to (Opi, a relativistic 
one selects a resonant mode of refractive index close to 
1, which corresponds to a frequency slightly larger than 
cOci- Thus whereas hybrid waves are strongly excited in a 
non-relativistic collisionless shock foot, lower extraordi- 
nary modes are moderately excited in a relativistic shock. 
They can nevertheless develop more intensively in case 
A and their instability corresponds to the micro-scale 
version with resonance of the previous compressive in- 
stability. Like for its MHD scale version, mostly heat- 
ing is expected from this instability; by heating we mean 
mixing of the streams in the foot region. 



OBSTRUCTION BY THE MEAN FIELD 

In this section we analyze the obstruction to the instabil- 
ity growth by the mean field (case B). Let us see first how 
it prevents Weibel instability. 

The condition that the pervading cosmic rays are 
not magnetized at the time scale the Weibel instability 
growth is naturally fulfilled since it requires 
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(23) 



A necessary condition for Weibel instability to devel- 
ope is that the ambient protons are non-magnetized at the 
time scale of the instability, which requires jisCOp* ^ (Olu, 
namely 
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(24) 



The most severe requirement is that the instability has 
time to grow significantly during the precursor crossing, 
namely: 

(25) 



which generally implies a very small mean field: 
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(26) 



When this condition is fulfilled, then the condition of 
non-magnetization of ambient protons is fulfilled. It has 
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FIGURE 4. Instability diagram: X = ^^r„ 

Y = r^S^sin^ OglATtncrnipV^ . Weibel instability develops for 

Y < 1, the whistler instability for < 1 and Y > 



been shown fie*] that the ambient electrons can be mag- 
netized; this does not quench Weibel instability. 

As for the whistler instability the analogous require 
ment is 
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Since k^5l can reach the value of the mass ratio mp/me, 
for Fj < mp/nie the whistler modes are less inhibited by 
the mean field than the Weibel waves. Nevertheless even 
with whistler mode excitation, the intensity of the mean 
magnetic field needs to be quite low in order that waves 
have time to grow, except, of course when the shock 
Lorentz factor is weak. 

Short scale magnetic fluctuations can undergo a fast 
growth in an ultrarelativistic shock precursor only if the 
ambient mean field is very weak; then the relativistic 
Fermi process is operative. The criteria is the following. 
Let X = ^F, and Y = I^Blsin^ eB/AKticrmpV^; then 
Weibel waves grow if F < 1; whistler modes grow if 
XY < 1 and the previous growth rate is valid for Y > 
X^ (otherwise there is no resonance); all the instability 
regimes are sketched in fig. 4. 

This is only in case A that we can expect a significant 
excitation of extraordinary ionic modes and magneto- 
sonic modes. They generate a pre-heating process in the 
precursor that concerns both electrons and protons. We 
cannot say at this stage of the investigation whether the 
incoming plasma is strongly heated in the shock foot. 
If electrons and protons are heated at similar relativistic 
temperature, it changes the instability dynamics. More- 
over it would start to slow down significantly the up- 
stream flow. 



Hot plasma in shock foot 

We do not know yet the details of a relativistic shock 
front. In the previous section we assumed that the shock 
front is structured like a non-relativistic front and just 
extended the non-relativistic results. Since MHD com- 
pressive instability and extraordinary ionic modes can be 
excited, we cannot exclude that the foot be full of rela- 
tivistically hot protons and electrons of similar temper- 
ature {frripC^). In that case the plasma response is dif- 
ferent, because the intermediate whistler range (and also 
extraordinary range) disappears; the plasma behaves like 
a relativistic pair plasma. In that case, only Weibel in- 
stability is expected and the limit on the intensity of the 
ambient magnetic field for the instability growth is such 
that 
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If this happens, this would be the ultimate criteria for 
developing Fermi acceleration at relativistic shocks. 



WHAT HAPPENS DOWNSTREAM? 

In the downstream plasma, the magnetic fluctuations 
generated by Weibel instability disappear rapidly be- 
cause they do not correspond to plasma modes. However 
whistler waves are transmitted and although they are not 
excited downstream, their damping is weak. When Fermi 
cycles develop, they create "inverted" distribution down- 
stream, that should produce a maser effect. Otherwise the 
intrinsic magnetic fluctuations would decay, only extrin- 
sic MHD turbulence would scatter large Larmor radius 
particles. 

Tangled magnetic field carried by the upstream flow 
are very compressed downstream and thus opposite 
polarization field lines come close together. This pro- 
duce magnetic reconnections in an usual regime where 
protons and electrons have a similar relativistic mass of 
order FsMpC^. Such a regime of reconnection deserves a 
specific investigation with appropriate numerical simula- 
tions. Despite magnetic dissipation, reconnections would 
probably create a chaotic flow that favors diffusion of 
particles from downstream to upstream. 



DISCUSSION 

Because the ambient mean magnetic field has a domi- 
nant transverse component at the relativistic shock front, 
the scattering of particles that are expected to undergo 
Fermi cycles requires very intense magnetic fluctuations 
at short scales. Short scales instabilities can be excited by 



the Fermi cycles, and eventually by the first Fermi cycle. 
However the mean transverse field again severely limits 
the efficiency of the instabilities by reducing drastically 
the length of the precursor. We found two regimes of effi- 
cient excitation of magnetic fluctuations suitable for scat- 
tering: one dominated by Weibel instability, the other by 
the resonant instability of the whistler modes. If particles 
are not too much pre-heated in the precursor, the latter is 
less inhibited by the mean field than the former. However 
it is also severely limited and no Fermi acceleration can 
occur, unless the mean field intensity is such that: 
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Now if electrons and protons are heated to relativistic 
temperature in the shock foot, only Weibel instability can 
make Fermi process to work provided that the mean field 
is so weak that 
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Indeed we found instabilities that can heat the precursor. 
As can be seen, the angle of the magnetic field with 
respect to the shock normal is a sensitive parameter. It 
turns out that the perpendicular collisionless shocks are 
unsteady in non-relativistic regime. They are probably 
unsteady in relativistic regime also. This is related to the 
development of whistler waves (f^S*] (2^). The shock 
reformation occurs at a time scale comparable to the 
particle crossing time in the foot. Such investigation for a 
relativistic shock would be important and would deserve 
heavy numerical simulations. 

Regarding the generation of UHECRs by relativistic 
shocks, the conclusions of this study is fairly pessimistic, 
because even if the Fermi cycles work with intense short 
scale magnetic fluctuations, the scattering time becomes 
longer and longer with the square of the particle energy. 
The Hillas criteria is no longer relevant; it should be re- 
placed by E„u,x ~ TecB{£cRy^^, where B is the quadratic 
average of the fluctuating magnetic field at the short scale 
ic- This new criteria is of no practical use, however it in- 
dicates that the very short scale £c makes the UHECR 
generation hopeless through this process. 

These topics suggest the development of heavy 3D 
PIC simulations for the nonlinear stage, as already 
started successfully ((HI, IH, [20]). 

This talk reveals the interest of developing commu- 
nications between three disciplines: High Energy Astro- 
physics, Particle Physics and Space plasma physics. 
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